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Figure  Captions 


1.  Simplified  tectonic  map  of  Eurasia.  Shield  locations  are  taken  from 
Burke  and  Sengdr  (1986)  and  the  Tethysides  complex  is  simplified  from 
Sengdr  (1984). 

2.  Map  of  events  used  in  this  study. 

3.  Map  of  stations  that  provided  data  for  this  study. 

4.  Comparison  of  seismograms  recorded  for  paths  through  high-Q  (top) 
and  low-Q  (bottom)  regions.  The  upper  seismogram  is  from  a 
magnitude  4.5  earthquake  in  the  Norwegian  Sea  (69.418N,  14.201E)  at 
a  distance  of  1 144  km  from  station  KONO  and  the  lower  seismogram  is 
from  a  magnitude  4.2  earthquake  in  Sichuan  Province,  China 
(29.996N,  99.338E)  at  a  distance  of  644  km  from  station  KMI. 

5.  Comparison  of  SSR  plots  obtained  from  the  seismograms  in  Figure  4. 
Note  that  the  value  obtained  for  frequency  dependence  (determined  from 
the  slope  of  the  least-squares  fit  line)  will  depend  somewhat  on  the 
selected  frequency  range. 

6.  Sampling  patterns  for  Lg  coda  used  in  this  study.  Each  Lg  coda  is 
assumed  to  sample  an  elliptical  area  corresponding  to  the  maximum 
lapse  time  used  in  the  analysis.  When  a  continental  boundary  is 
encountered  we  assume  it  is  a  barrier  to  Lg  propagation  and  therefore 
forms  part  of  the  boundary  for  the  sampling  area. 

7.  Tomographic  map  of  Lg  coda  Q  at  1  Hz  for  most  of  Eurasia. 

8.  Tomographic  map  of  the  frequency  dependence,  Tj,  of  Lg  coda  Q  at  1 
Hz  for  most  of  Eurasia. 

9.  Point  spreading  function  (psf)  plots  for  grids  centered  at  70.93ON  - 
41.960E,  at  63.790E  -  125.530E,  at  48.770N  -  27.850E,  and  at 
32.40ON  -  96.830E. 

10.  Error  distribution  for  the  Qo  image  of  Eurasia.  The  error  is  tested 
empirically,  based  upon  the  sample  standard  error  in  the  Qo  and  Tj 
values  calculated  from  real  data  and  the  pseudo-random  sign  generators. 
See  text  for  details. 

1 1 .  Error  distribution  for  the  T|  image  of  Eurasia.  The  error  is  tested 
empirically  as  described  in  the  caption  for  Figure  10. 
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Introduction 

Lg  is  often  the  most  prominent  and  robust  phase  on  short-period  and  broad-band 
seismograms  recorded  in  continental  regions.  For  that  reason  it  has  been  used  to  determine 
magnitudes  of  small  events  at  regional  distances  (Baker,  1970;  Nuttli,  1973),  yields  of 
nuclear  explosions  (Nuttli,  1986),  and  regional  variations  of  crustal  attenuative  properties 
(Campillb  &  Plantee,  1991;  Mitchell,  1981;  Sutton,  Mitronovas,  &  Pomeroy,  1967).  It  can 
be  treated  as  a  superposition  of  higher-mode  Rayleigh  waves  that  propagate  in  the  continent 
crust  (Kennett,  1984;  Knopoff,  Schwab,  &  Kausel,  1973;  Panza  &  Calcagnile,  1975).  The 
onset  of  Lg  arrives  with  a  group  velocity  of  3.3  -  3.6  km/s,  with  young  active  regions 
being  slower  than  stable  shields. 

A  coda  of  up  to  several  minutes  duration  follows  Lg,  the  predominant  frequency  of 
which  is  lower  than  the  direct  phase,  and  decreases  with  time.  The  arrival  of  Lg  coda 
energy  is  omnidirectional,  indicating  that  it  is  scattered  from  heterogeneities  in  the  crust 
(Der,  Marshall,  O'Donnell,  &  McElfresh,  1984).  Being  a  scattered  wave,  it  is  sensitive  to 
properties  in  a  volume  of  the  crust  that  surrounds  the  great  circle  path,  including  both  the 
seismic  source  and  the  receiver,  rather  than  to  properties  along  a  single  ray  path  between 
them.  Scattering  had  earlier  been  inferred  as  the  cause  for  coda  following  body  waves  at 
short  distances  (Aki,  1969;  Aki  &  Chouet,  1975). 

The  attenuation  of  both  Lg  and  Lg  coda  is  often  expressed  as  the  inverse  of  a  quality 
factor,  Q.  In  the  present  paper,  we  use  Q^g  and  Q^g  to  denote,  respectively,  the  quality 

factors  of  the  direct  Lg  phase  and  Lg  coda.  An  advantage  that  use  of  Lg  coda  has  over  use 
of  the  direct  Lg  phase  for  studying  regional  variations  of  attenuation  is  that  measurements 
can  be  made  using  data  from  only  a  single  station  and  do  not  require  knowledge  of  the 
source  spectrum.  Q  determinations  from  direct  Lg  require  either  a  single-station  recording 
from  an  event  with  a  known  source  spectrum  (or  a  spectrum  that  can  be  determined)  or 
records  from  two  or  more  stations  that  lie  on  a  great-circle  path  with  a  seismic  event.  Q^g 
determinations  are,  therefore,  usually  much  more  plentiful  than  those  of  Qlj  and  provide 
better  areal  coverage  in  studies  of  attenuative  properties. 

A  tomographic  method  for  studying  regional  variations  of  Lg  coda  Q  was  developed  by 
(Xie  &  Mitchell,  1990a)  and  was  applied  to  Africa  for  frequencies  near  1  Hz.  It  assumes 
that  the  area  sampled  by  Lg  coda  forms  an  ellipse  with  the  source  at  one  focus  and  the 
recording  station  at  the  other  (Malin,  1978).  The  areal  coverage  of  the  ellipse  for  a  given 
source-station  pair  increases  with  increasing  lapse  time.  Xie  and  Mitchell  (1990a)  applied  a 
back-projection  inversion  technique  (Dines  &  Lytle,  1979)  to  obtain  maps  of  the  variation 
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of  Qo  and  r\  ( Q^g  at  1  Hz  and  its  frequency  dependence)  for  all  of  Africa  from  many  event- 

station  pairs.  We  present  details  of  the  method  in  later  sections. 

The  African  study  revealed  striking  variations  in  Q^g,  with  low  values  characterizing 

regions  in  which  tectonic  activity  has  occurred  since  the  Mesozoic  era.  These  include  the 
East  African  rift  zone  (the  most  prominent  feature  on  the  Q  map),  the  Atlas  Mountains,  the 
Cape  fold  belt,  and  the  Cameroon  line.  High  Q^g  values  characterize  the  central  portions  of 

the  cratons  which  came  together  to  form  Africa  about  550  M.y.  ago.  Mitchell  (1995) 
recently  suggested  that  intermediate  values  of  Q^g  that  occur  in  the  fold  belts  between  the 

cratons  reflect  tectonic  activity  that  occurred  during  the  formation  of  the  African  continent  in 
the  early  Paleozoic  era.  A  revised  map  of  frequency  dependence  values  for  Africa  (Mitchell, 
1995)  indicates  that  T|  varies  inversely  with  Q^g,  with  low  frequency  dependence  being 

characteristic  of  the  cratons  and  high  frequency  occurring  in  the  low-Q  regions. 

This  report  applies  the  methods  of  Xie  and  Nuttli  (Xie  &  Nuttli,  1988)  and  (Xie  & 
Mitchell,  1990a)  to  Eurasia,  a  composite  continent  consisting  of  numerous  diverse  blocks 
which  began  a  complex  coalescence  during  the  Devonian  period,  about  370  million  years 
ago  (Zonenshain,  Kuzmin,  &  Natapov,  1990).  In  contrast  with  Africa,  for  which  the 
tectonics  is  dominated  by  two  events  (the  accretion  of  the  cratons  during  the  Pan- African 
tectono-thermal  event  about  550  M.y.  ago  and  the  recent  extension  in  the  East  African  rift), 
Eurasia  has  undergone  numerous  and  diverse  periods  of  tectonic  and  organic  activity. 
Moreover,  whereas  Africa  is  characterized  by  relatively  low  relief  and  modest  sedimentary 
accumulations,  Eurasia  contains  both  the  world’s  highest  and  lowest  points  and  several  of 
its  greatest  sedimentary  basins.  Later  sections  present  maps  of  Lg  coda  Q  and  its  frequency 
dependence  at  frequencies  near  1  Hz  for  most  of  Eurasia.  We  then  compare  them  with 
results  from  more  localized  studies  of  QLg  and  Q^g  as  well  as  with  fundamental-mode 

surface  attenuation,  and  correlate  them  with  tectonic  features  of  Eurasia. 

Tectonic  Setting  and  Evolution  of  Eurasia 

Eurasia  is  a  collage  of  continental  blocks  that  have  merged  over  a  long  period  of  time 
beginning  in  the  early  Paleozoic  era.  The  two  shields  of  northern  Eurasia,  the  East 
European  Shield  and  the  Siberian  Shield,  are  believed  to  have  amalgamated  from  smaller 
continental  blocks  and  to  have  become  fully  formed  by  about  1700  Ma  (Zonenshain,  et  al., 
1990).  The  East  European  block  collided  with  North  America  in  the  early  Paleozoic  to  form 
Laurentia,  an  event  recorded  by  the  Caledonian  fold  belt  in  Scandinavia,  Scotland,  and 
Svalbard.  Laurentia  and  the  Siberian  block  remained  separated  from  one  another  and  from 
numerous  smaller  cmstal  blocks  until  convergence  began  in  the  late  Devonian  period. 
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Zonenshain  et  al.  (1990)  discuss  the  formation  of  Eurasia,  with  emphasis  on  the 
evolution  of  the  crust  within  the  former  USSR,  while  Sengor  (1984;  1987)  details  the 
evolution  of  the  Tethy sides,  the  orogenic  complex  that  extends  from  Spain  to  eastern  Asia 
forming  the  Alpine-Himalayan-Indonesian  mountain  ranges.  Many  major  features  related  to 
the  evolution  of  Eurasia  appear  in  the  simplified  tectonic  map  in  Figure  1.  Names  of 
tectonic  features  on  that  map  are  mainly  taken  from  Sengor  (1984)  and  Burke  (1986)  but 
we  replace  their  terms  Serindia  and  Angara  Shield,  respectively,  with  the  more  widely 
known  terms  Tarim  Basin  and  Siberian  Shield. 

The  East  European  and  Siberian  shields  collided  with  the  Altaid  complex  and  with  each 
other  beginning  about  310  Ma.  Accretion  of  numerous  smaller  plates  culminated  in  a  single 
Laurasia  about  280  Ma.  The  Altaids  were  cemented  between  the  Siberian  and  East 
European  plates  and  orogeny  occurred  at  that  time  forming  the  Urals  at  the  western  margin 
of  the  Altaids.  At  this  time  convergence  also  began  between  the  Paleo-Tethys  oceanic  plate 
and  Laurasia  along  a  6000  km  long  zone.  As  the  Paleo-Tethys  was  consumed  beneath  the 
Eurasian  margin,  parts  of  the  newly  formed  Neo-Tethys  (separated  from  the  Paleo-Tethys 
by  the  Cimmerian  continent)  widened  over  some  of  its  length. 

At  220  Ma  the  plates  were  extensively  reorganized,  with  the  collision  of  North  China 
block  with  Eurasia  and  the  nearly  complete  closure  of  Paleo-Tethys.  Compression  ceased 
within  Eurasia  and  extension  affected  parts  of  western  Siberia  and  the  Barents  Sea  Basin. 

The  margin  of  the  Africa- Arabia  plate  began  to  progressively  converge  with  Eurasia  at 
about  160  Ma  gradually  closing  the  Tethy s  ocean.  The  breakup  of  Pangea  occurred 
between  160  and  130  Ma  while  Africa- Arabia  was  converging  with  Laurasia.  Eurasia 
separated  from  North  America  at  about  80  Ma  while,  at  the  same  time,  exotic  terrains  from 
the  Pacific  continued  to  be  deposited  at  the  eastern  margin  of  the  plate. 

During  the  Cenozoic  (after  65  Ma)  plate  convergence  began  between  Eurasia  and  India 
in  the  south,  and  continued  between  Eurasia  and  Africa-Arabia  in  the  southwest,  and  the 
Pacific  plate  in  the  east.  Simultaneously,  the  width  of  the  Tethys  substantially  decreased. 

The  complex  evolutionary  history  of  Eurasia  has  left  a  geologically  diverse  continent 
with  much  lateral  complexity.  In  addition  to  features  formed  at  the  margins  of  the  colliding 
plates  during  the  long  evolutionary  history  of  Eurasia,  some  cratonic  regions  have  been 
covered  by  younger  sediments  that  may  reach  great  thicknesses.  Examples  of  deep 
sedimentary  basins  include  those  in  the  Barents  shelf  (Levshin  &  Berteussen,  1979)  where 
sediment  depths  reach  about  25  km  and  the  Peri-Caspian  basin  (Neprochnov,  1968)  where 
they  are  as  thick  as  20  km. 
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Figure 


Data  Coverage 

Data  for  this  study  were  generated  by  earthquakes  located  within  or  just  at  the  boundaries 
of  the  Eurasian  plate.  We  avoided  earthquakes  outside  those  boundaries  because  of  the 
likelihood  of  systematic  errors  in  Q  and  its  frequency  dependence  known  to  be  produced 
when  the  scattering  ellipse  associated  with  any  event-station  pair  includes  a  significant 
portion  of  oceanic  crust  (Xie  &  Mitchell,  1990b).  Because  of  the  heterogeneity  of 
earthquake  occurrence  throughout  Eurasia,  there  were  many  more  usable  events  in  the 
southern  portion,  especially  in  China  and  the  Middle  East  than  in  the  northern  portion  of  the 
continent  (Figure  2). 

Recording  stations  used  in  this  study  were  widely,  and  fairly  evenly,  distributed  across 
Eurasia  (Figure  3).  Data  were  obtained  from  stations  belonging  to  seven  different  networks 
of  stations  (Table  1),  some  of  which  provided  digital  data  and  some  of  which  provided  data 
in  analog  form.  Digital  data  were  available  from  networks  operated  by  Incorporated 
Institutions  for  Seismology  (IRIS),  the  Chinese  Digital  Seismic  Network  (CDSN),  the 
Seismic  Research  Observatories  (SRO),  the  Digital  World-Wide  Standard  Seismograph 
Network  (DWWSSN),  and  the  Norwegian  Seismic  Array  (NORSAR),  while  data  from 
stations  in  the  former  Soviet  Union  (MOS),  and  from  the  World-Wide  Standard 
Seismograph  Network  (WWSSN)  were  in  analog  form. 

All  of  the  Lg  coda  Q  determinations  reported  in  this  study  are  newly  determined  except 
those  for  the  Indian  subcontinent  (John,  1983)  and  the  Arabian  peninsula  (Ghalib,  1992). 
which  were  obtained  in  earlier  studies.  We  analyzed  a  few  traces  in  the  Arabian  Peninsula 
using  the  SSR  method  and  found  them  to  be  consistent  with  the  earlier  obtained  values. 

This  configuration  of  events  and  stations  in  Figures  2  and  3  produced  375  seismograms 
from  298  seismic  events  and  provided  coverage  for  all  of  Eurasia  except  for  Spain, 
Portugal,  Italy  and  Ireland  in  Europe  and  for  portions  of  north  easternmost  and  south 
easternmost  Asia  and  the  southern  tip  of  India.  The  regions  in  Europe,  southeastern  Asia, 
and  India  lack  coverage  because  it  is  not  possible  obtain  data  there  that  is  not  affected  by  the 
nearby  oceans,  whereas  the  gap  in  northeastern  Siberia  occurs  because  of  a  lack  of  station 
coverage.  Data  coverage  across  the  rest  of  Eurasia  is,  of  course,  not  uniform  because  of  the 
uneven  distribution  of  events.  The  best  coverage  is  available  for  a  band  that  stretches  from 
China  to  the  Middle  East  because  of  the  large  number  of  events  and  the  relatively  good 
station  coverage  there. 

The  character  of  Lg  and  its  coda  varies  greatly  for  different  paths  in  Eurasia.  Figure  4 
displays  two  seismograms,  one  from  station  KONO  in  southern  Norway  produced  by  an 
earthquake  in  the  Norwegian  Sea  and  the  other  from  station  KMI  in  southern  China 
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Figure  4 


Table  1.  Seismographic  Stations  Used  in  This  Study 


Station 

Latitude 

Longitude 

Network 

ABKT 

37.9304°N 

58.1189°E 

IRIS 

AKU 

56.4302°N 

58.5625°E 

IRIS 

GAR 

39.0000°N 

70.3167°E 

IRIS 

GNI 

40.0530°N 

58.1189°E 

IRIS 

KBS 

78.9256°N 

11.9417°E 

IRIS 

KTV 

43.9553°N 

42.6863°E 

IRIS 

KONO 

59.6491°N 

9.5982°E 

IRIS 

NRIL 

69.5049°N 

88.4414°E 

IRIS 

NVS 

54.8404°N 

83.2346°E 

IRIS 

OBN 

55.1138°N 

36.5687°E 

IRIS 

TLY 

51.6807°N 

103.6438°E 

IRIS 

BJI 

40.0403'='N 

116.1750°E 

CDSN 

HIA 

49.2667°N 

119.7417°E 

CDSN 

KMI 

25.1233°N 

102.7400°E 

CDSN 

LZH 

36.0867°N 

103.8444°E 

CDSN 

MDJ 

44.6164°N 

129.5919°E 

CDSN 

WMQ 

43.8211°N 

87.6950°E 

CDSN 

ANTO 

39.8689°N 

32.7936°E 

SRO-sp 

CHTO 

18.8138°N 

98.9438°E 

SRO-sp 

GRFO 

49.6919'“N 

11.2217°E 

SRO-sp 

SHIO 

25.5666°N 

91.8830“E 

SRO-sp 

APA 

67.5500°N 

33.3333°E 

MOS 

BOD 

57.8499'>N 

114.1833°E 

MOS 

KHE 

80.6167°N 

58.0500'='E 

MOS 

TIK 

71.6333°N 

128.8667°E 

MOS 

TUP 

54.4333°N 

119.9000°E 

MOS 

UZH 

48.6333°N 

22.3000°E 

MOS 

EIL 

29.5500°N 

34.9500°E 

WWSSN-sp 

JER 

31.7719°N 

35.1972°E 

WWSSN-sp 

NDI 

28.6833°N 

77.2167°E 

WWSSN-sp 

NIL 

33.6500°N 

73.2517°E 

WWSSN-sp 

NUR 

60.5090°N 

24.6514°E 

WWSSN-sp 

POO 

18.5333°N 

73.8500°E 

WWSSN-sp 

QUE 

30.1883°N 

66.9500°E 

WWSSN-sp 

SHI 

29.6445°N 

52.5261°E 

WWSSN-sp 

SHL 

25.5667°N 

91.8833°E 

WWSSN-sp 

TAB 

38.0675°N 

46.3267°E 

WWSSN-sp 

UME 

63.8150°N 

20.2367°E 

WWSSN-sp 

BER 

60.3869°N 

5.3258°E 

DWWSSN-sp 

KEY 

69.7553°N 

27.0067°E 

DWWSSN-sp 

06C 

60.7473°N 

11.4583°E 

NORSAR 
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produced  by  an  earthquake  in  western  China.  The  Lg  onset  is  roughly  at  2  minutes  from 
the  plot  origin  in  both  cases.  The  predominant  frequency  is  much  higher  for  the  KONO 
record  than  for  the  KMI  record,  even  though  the  event-station  distance  between  the 
recorded  earthquake  and  KONO  is  much  longer  (1 144  km)  than  the  event-station  distance 
between  the  recorded  earthquake  and  KMI  (644  km).  This  obvious  depletion  of  high 
frequencies  at  KMI  relative  to  that  at  KONO  indicates  that  energy  is  attenuated  much  more 
severely  in  China  than  it  is  in  Norway.  The  large  difference  in  frequency  content  suggests, 
moreover,  that  variations  in  the  attenuative  properties  in  different  regions  of  Eurasia  can  be 
easily  mapped. 

The  development  of  Lg  coda  relative  to  that  of  direct  Lg  for  the  seismograms  of  our 
study  is  highly  variable.  Direct  Lg  is  often  large  and  sharp,  in  which  case  it  is  always 
followed  by  a  coda,  but  other  times  it  is  very  emergent,  sometimes  not  easily  observed  in 
noise.  In  that  latter  case  usable  Lg  coda  is  often  present,  but  has  an  emergent  beginning  and 
a  gradual  die-off  with  time 

Data  Processing 

Xie  and  Nuttli  (1988)  developed  a  stochastic  model  for  the  interpretation  of  single-trace  Lg 
coda  waves  and  an  inversion  method,  based  on  that  model,  to  determine  Lg  coda  Q  and  its 
frequency  dependence  at  1  Hz  (Qo  and  T|,  respectively).  Advantages  of  their  inversion 
method  are  that,  with  data  of  sufficient  quality,  tradeoffs  between  Qo  and  T|  are  avoided 
(often  a  problem  with  earlier  methods)  and  that  variances  in  those  quantities  are  reduced 
compared  to  other  inversion  methods.  Application  of  their  stochastic  model  to  examples  of 
recorded  ground  motion  showed  that  the  randomness  that  occurs  in  Lg  coda  waves  can  be 
approximated  by  simple  band-limited  white  Gaussian  noise,  a  result  that  allows  us  to  make 
quantitative  estimates  of  variances  when  estimating  Lg  coda  attenuation. 

The  inversion  method,  termed  the  stacked  spectral  ratio  (SSR)  method,  has  been 
discussed  elsewhere  (Xie  &  Mitchell,  1990a;  Xie  &  Nuttli,  1988);  for  that  reason  it  is  only 
summarized  here.  It  begins  by  dividing  the  Lg  coda  time  series  into  a  number,  Nw,  of  non¬ 
overlapping  time  windows  of  constant  length  and  applying  the  discrete  Fourier  transform  to 
each  window  to  obtain  the  geometrical  mean  of  the  amplitude  spectrum  at  each  frequency. 
The  length  of  each  group  velocity  window  for  Lg  is  given  by 

C/„  =  VT  (1  /  V  •  —  1  /  )  (1) 

where  Tm  is  the  mean  lag  time  for  the  mth  window,  and  v,  Vmax.  and  Vmin  are, 
respectively,  the  average,  maximum,  and  minimum  group  velocities  of  the  coda  wave  train. 
Geometrical  spreading  is  obtained  using 

G^={2KRr'\v^iil  (2) 
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where  R  is  the  event-station  distance.  The  SSR  is  defined  as 

M  ^ 


1  M 


^  ^m) 


log 


Gffl  V^M+ra  {-At  ,m+Af  ^ 


(3) 


where  M  is  Nw/2  when  Nw  is  even  and  (Nw  +  l)/2  when  Nw  is  odd.  All  of  the  quantities 
on  the  right-hand  side  of  (3)  are  calculable  from  a  coda  time  series. 

Assuming  that  the  frequency  dependence  of  Lg  coda  Q  can  be  expressed  by  a  power 
law,  Q(f)  =  Qof",  where  Qo  and  r\  are,  respectively,  Q  and  its  frequency  dependence  at  1 

Hz,  then 

logF^  =  (1  -  7/)log/^  - log(2„  +  e  (4) 

where  e  represents  a  random  error.  Theoretical  and  observational  calculations  (Xie  & 
Nuttli,  1988)  have  shown  that  the  SSR  provides  a  statistically  stable  Q  estimator,  with 
the  standard  error  in  Qo  usually  being  an  order  of  magnitude  smaller  than  Qo  itself. 

The  SSR  method  was  applied  to  all  traces  beginning  at  a  group  velocity  of  3.15  km/s 
and  continuing  on  in  time  to  cover  a  window  length  of  300  seconds  or  more,  with  the 
average  length  being  about  450  seconds  and  the  maximum  duration  being  600  seconds.  We 
tried  to  minimize  variations  in  coda  length  from  path  to  path  in  order  to  avoid  possible 
systematic  biases  that  might  be  proportional  to  coda  duration.  Using  different  window 
lengths  for  the  same  seismogram  we  found  that  Lg  coda  Q  sometimes  increased  as  much  as 
20%  if  coda  length  was  doubled.  A  similar  dependence  on  coda  length  was  noted  in  an 
earlier  study  using  shear  waves  at  shorter  distances  (Pulli,  1984). 

Another  possible  source  of  bias  in  our  determinations  can  arise  if  a  large  portion  of  the 
scattering  ellipse  overlaps  oceanic  regions  (Xie  &  Mitchell,  1990b).  For  that  reason,  we 
avoided  events  that  occurred  in  oceanic  regions  as  well  as  continental  events  for  which  a 
large  fraction  of  the  scattering  ellipse  overlaps  an  oceanic  region. 

Figure  5  shows  example  SSR  plots  for  the  seismograms  in  Figure  4.  The  plots  indicate 
that  for  station  KONO  Qo  (Qlg^t  1  Hz)  is  770±274  and  the  frequency  dependence  (from 

the  slope  of  the  least-squares  fit)  is  0.12±0.20,  whereas  for  station  KMI  Qo  is  234±66  and 
frequency  dependence  is  0.17±0.01.  The  plots  indicate  that  if  SSR  values  are  linear  with 
frequency  over  a  significant  range  around  the  frequency  of  interest,  then  reliable  values  of 
both  Qo  and  r\  can  be  obtained. 


Mapping  Lg  Coda  Q  Variations  -  Method 

Although  seismic  tomography  has  been  very  successful  in  mapping  velocities  in  the 
Earth,  it  has  not  been  used  extensively  to  map  variations  of  attenuative  properties.  The  first 
such  study  for  regional  phases  was  done  for  the  African  continent  where  a  set  of  Lg  coda 
records  that  covered  much  of  the  continent  was  assembled  and  the  SSR  method  was 
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applied  to  them  (Xie  &  Mitchell,  1990a).  The  resulting  tomographic  image  revealed  striking 
regional  variations  of  Lg  coda  Q  that  coincide  with  some  of  the  major  tectonic  features  of 
Africa. 

In  the  present  study  we  apply  the  same  method  to  Eurasia.  It  is  a  back-projection 
algorithm  that  produces  a  tomographic  image  of  Q^g  over  any  broad  region  using  a  number 

(Nc)  of  individual  determinations  of  Q[g  from  Lg  coda  recordings  within  that  region. 


Since  Lg  coda  consists  of  scattered  waves  the  energy  arriving  at  a  recording  station  does 
not  follow  a  straight  path  between  from  the  seismic  event,  but  is  contained  within  a  volume 
of  the  crust  that  includes  the  source  and  recording  station. 

For  single  scattering  the  surface  projection  of  the  scattered  energy  takes  the  form  of  an 
ellipse  with  the  source  at  one  focus  and  the  recording  station  at  the  other  (Malin,  1978). 
The  ellipses  cover  larger  and  larger  areas  as  lag  time  of  the  coda  increases.  The  elliptical 
areas  for  the  event-station  pairs  used  in  Eurasia  are  shown  in  Figure  6  at  the  maximum 
lapse  time  used  for  the  SSR  analysis  described  earlier.  If  coda  waves  are  scattered  through 
a  volume  as  approximated  by  the  elliptical  surface  projections  in  Figure  6,  and  if  many 
events  and  recording  stations  are  available,  then  much  overlap  occurs  among  some  of  the 
ellipses.  This  overlap  provides  data  redundancy  that  is  ideal  for  applying  modem  inversion 
techniques  to  map  regional  variations  of  attenuation.  The  process  for  doing  that  is  described 
in  the  following  paragraphs. 

To  apply  the  tomographic  method  we  divide  the  region  of  interest  (in  the  present  case, 
most  of  Eurasia)  into  a  number,  Nc,  of  cells  that  are  small  enough  to  allow  us  to  resolve 
large-scale  tectonic  features.  Qm,  Lg  coda  Q  within  the  mth  cell,  is  assumed  to  be  constant 
throughout  the  cell.  For  Eurasia,  we  take  the  cell  dimensions  to  be  3°  x  3°.  Qn,  the  Lg  coda 
Q  value  calculated  for  the  nth  time  series,  gives  the  areal  average  of  Lg  coda  Q  over  the 
elliptical  area  sampled  by  the  coda  waves  received  at  the  maximum  lapse  of  the  Lg  coda 
waves  of  interest.  If  we  denote  the  area  that  the  nth  ellipse  overlaps  the  mth  grid  by  Smn, 
then 


Qn 


Qn 


n=l,  2,  .  .  .,Nd 


(5) 


where 


\  =  S-'*  (6) 

7=1 

and  8n  is  the  residual  due  to  the  errors  in  the  measurement  and  modeling  Lg  coda.  We  solve 
(5)  by  a  back-projection,  or  ART  technique  (Gordon,  1974),  that  has  been  used 
successfully  in  seismic  velocity  tomography  (Humphreys  &  Clayton,  1988;  McMechan, 
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1983;  Suetsugu  &  Nakanishi,  1985).  The  numerical  procedure  for  that  solution  is 
described  by  Xie  and  Mitchell  (1990a)  and  will  not  be  repeated  here. 

A  disadvantage  of  the  back-projection  technique  is  that  it  does  not  provide  formal 
estimates  of  resolution  or  standard  error.  Humphreys  and  Clayton  (1988)  have,  however, 
suggested  using  a  “point  spreading  function”  as  an  approximation  for  the  resolution  kernel. 
The  point  spreading  function  (hereafter  abbreviated  psf)  is  obtained  by  constructing  a 
model  in  which  has  the  value  of  unity  in  a  single  cell  at  the  geographical  point  of 
interest  and  is  zero  everywhere  else.  We  then  compute  synthetic  data  using  this  model  and 
perform  an  inversion.  The  resulting  image  gives  the  psf  at  the  point  of  interest,  a  quantity 
whose  area  and  falloff  with  distance  gives  a  measure  of  resolution. 

The  effect  that  random  noise  in  our  Lg  coda  has  on  our  tomographic  images  can  be 
empirically  tested  using  the  sample  standard  error  in  Qn  caused  by  randomness  of  the 
SSR's  (Xie  &  Nuttli,  1988).  If  the  standard  error  for  each  determination  of  Qn  is  denoted 
by  6Qn,  we  can  construct  a  noise  series  for  each  Lg  coda  time  series  that  has  an  absolute 
value  equal  to  5Qn  and  a  sign  that  is  chosen  randomly.  The  nth  term  of  the  noise  series  is 
added  to  Qn  and  the  sum  of  the  two  series  are  then  inverted  to  obtain  a  new  Qm  image.  The 
difference  between  this  image  and  the  original  one  gives  us  an  error  estimate  for  the  Qm 
values.  Since  the  sign  of  5Qn  is  chosen  by  a  random  binary  generator,  the  process  should 
be  repeated  several  times  to  obtain  an  average  error  estimate. 

Mapping  Lg  Coda  Q  Variations  -  Results 

We  used  the  back  projection  method  described  above  to  image  lateral  variations  of  Lg  coda 
Q  across  most  of  Eurasia.  The  area  for  which  data  were  available  extends  between  about 
15°  N  and  nearly  85°  N  and  between  about  10°  W  and  about  160°  W.  The  map  depicting 
Lg  coda  scattering  ellipses  (Figure  6)  indicates  the  extent  of  our  data  coverage  for  the 
Eurasian  plate.  Northern  Eurasia  is  at  high  latitudes  where  the  area  of  any  n-degree  by  n- 
degree  cell  covers  a  much  smaller  than  it  does  at  lower  latitudes.  In  order  that  the  scattering 
ellipses  cover  approximately  the  same  number  of  cells  at  both  and  low  latitudes,  before 
dividing  Eurasia  into  cells,  we  performed  a  coordinate  transformation  to  center  Eurasia  at 
the  equator.  After  that  transformation,  we  divided  the  area  of  study  into  3°  by  3°  cells  and 
inverted  our  Qo  values  to  obtain  a  tomographic  image  of  Eurasia.  Following  the  inversion, 
we  transformed  the  image  to  its  original  position.  Although  this  process  produces 
approximately  uniform  weighting  for  all  cells  in  Eurasia,  it  distorts  the  image  at  high 
latitudes. 

Lateral  variations  of  Qo  and  r) 
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For  our  inversion  process,  using  equation  5,  we  used  a  starting  a  starting  model  in 
which  Q^=  600,  m  =  1,2,  ...  ,  Ng  for  the  entire  region  of  study.  The  equation  was 
inverted,  after  substituting  in  values  for  the  various  areas  and  determined  Qo  values  for  Qn. 
Residuals  were  calculated  for  each  iteration  and  the  stopping  criterion  was  either  that  the 
number  of  iterations  is  20  or  the  relative  change  in  the  sum  of  squares  of  residuals  is  less 
than  10‘5.  Figure  7  displays  the  tomographic  image  of  Qo  for  Eurasia.  Qo  ranges  between 
about  200  and  about  1000,  with  the  high  values  being  centered  at  three  widely  separated 
points. 

The  process  for  determining  the  frequency  dependence  of  Lg  coda  Q  begins  by  using 
the  Qo  and  T|  values  obtained  for  each  trace  to  estimate  Lg  coda  Q  at  another  frequency;  that 
being  3  Hz  in  this  study.  These  new  values  were  used  as  Qn  in  equation  5  and  a  map  Lg 
coda  Q  at  3  Hz  was  then  obtained.  Using  the  maps  of  Lg  coda  Q  at  1  Hz  and  3  Hz  we 
calculate  the  frequency  dependence  of  Lg  coda  Q  using  the  relationship 


Figure  8  shows  the  distribution  of  Tj  across  Eurasia  resulting  from  this  process. 

Figures  7  and  8  do  not  suggest  any  clear  relationship  between  Qo  and  Tj  for  Eurasia  as 
has  been  suggested  by  (Nuttli,  1988).  Possible  relationships  of  those  quantities  to  the 
tectonics  of  Eurasia  are,  however,  discussed  in  the  next  major  section. 

Resolution  and  Error 

Figure  9  presents  psf  plots  for  four  selected  regions  of  Eurasia.  As  discussed  earlier, 
the  areas  covered  by  these  plots  and  the  rapidity  of  their  falloff  with  distance  from  the 
center  cell  can  be  used  to  estimate  our  ability  to  resolve  features  of  the  mapped  Qo  and  Tj 
distributions.  Of  the  four  psf  plots  the  best  resolution  occurs  in  China  where  data  coverage 
in  densest  and  the  poorest  resolution  is  in  northeastern  Asia  where  coverage  is  the  poorest. 

The  method  described  in  the  previous  section  permits  estimates  of  the  effect  of  random 
errors  in  Lg  coda  Q  measurements  on  our  images  of  Qo  and  Tj.  Five  tests  were  run  to 
estimate  the  error  in  the  images  of  the  two  quantities.  In  each  of  the  tests  we  constructed 
two  noise  series,  in  both  of  which  the  absolute  values  of  the  nth  terms  equal  the  standard 
errors  found  for  Qo  and  Tj  from  the  nth  coda  seismogram.  The  signs  of  the  nth  terms  were 
randomly  generated  by  a  random  binary  generator.  The  nth  terms  of  the  two  noise  series 
were  then  added  to  the  nth  measurements  of  Qo  and  Tj  measurements  to  construct  new 
synthetic  maps  of  Qo  and  Tj.  The  differences  between  these  new  maps  and  the  original  ones 
were  calculated  and  stored  after  each  test.  The  average  values  of  the  differences  for  the  five 
tests  give  empirical  estimates  of  absolute  error  in  Qo  and  Tj  as  shown  in  Figures  10  and  1 1. 
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The  largest  errors  in  both  the  Q©  and  T)  maps  largely  coincide  with  the  higher  Qo  values 
in  Figure  10.  Fractional  errors  are,  however,  with  some  exceptions,  relatively  uniform 
throughout  Eurasia.  It  is  perhaps  noteworthy  that  the  errors  are  quite  high  in  a  region 
within  the  East  European  Shield,  where  both  Qo  and  ri  are  high,  a  situation  that  contrasts 
with  the  expected  relationship  between  Qo  and  rj  (Nuttli,  1988). 

Comparison  of  mapped  Qo  with  values  obtained  in  regional  studies 

Several  determinations  of  Lg  Q  and  Lg  coda  Q  have  been  made  in  localized  studies  in 
Eurasia.  In  places  where  both  quantities  have  been  determined  they  commonly  take  on 
similar  values  (Der,  et  al.,  1984;  Nuttli,  1988;  Xie  &  Nuttli,  1988).  For  that  reason  we 
have  compared  Qo  values  determined  in  the  present  study  with  regionally  derived  values  for 
both  Lg  Q  and  Lg  coda  Q.  The  Qo  and  Tj  values  in  eastern  China  (Figures  7  and  8)  are 
similar  to  those  obtained  in  a  localized  study  of  Lg  coda  Q  there  (Chun,  Liu,  Zhu,  &  Shih, 
1994).  Similarly,  regional-scale  studies  of  Lg  Q  in  the  Barents  shelf  (Nuttli,  1988),  in  Iran 
(Nuttli,  1980),  central  Asia  (Xie,  Cong,  &  Mitchell,  1996),  and  the  Tibetan  Plateau 
(McNamara,  Owens,  &  Walter,  1996)  have  yielded  values  that  are  consistent  with  those  of 
the  continent-scale  measurements  (Figure  7). 

Qo  and  Tectonics 

Figure  7  shows  three  regions  where  Qo  is  higher  than  in  surrounding  regions.  One  is  a 
broad  band  that  extends  from  Scandinavia  and  the  Barents  shelf  through  all  of  northern 
Asia  except  for  one  roughly  circular  region  centered  to  the  southeast  of  Novaya  Zemlya. 
This  band  coincides  with  the  two  shields  of  northern  Eurasia,  the  Ural  Mountains  and  the 
Altaid  belt  between  them.  A  long  broad  band  with  low  Qo  extends  from  the  southern  part 
of  western  Europe,  through  the  Middle  East  and  central  Asia,  into  China  and  southeastern 
Siberia. 

Most  of  the  low  Q's  in  Figure  7  correspond  with  the  Tethysides  tectonic  region  of 
Figure  1,  but  they  also  cover  the  Arabian  Peninsula,  a  shield  region,  where  high  Q's  are 
normally  expected.  The  Tethysides  is  a  deformed  region  that  began  forming  about  280  My 
ago  when  the  Paleo-Tethys  Ocean  began  to  be  consumed  beneath  Eurasia.  The  Alpinotype 
Tethysides  (Figure  1)  form  the  southern  portion  of  the  Tethysides  and  have  been  deformed 
as  a  consequence  of  orogeny  while  the  Germanotype  Tethysides  are  a  more  northerly  band 
that  consists  primarily  of  dominantly  brittle,  blocky  structures  such  as  grabens,  block 
uplifts,  and  large  strike-slip  fault  systems  (Sengor,  1984). 

Another  shield  where  our  map  exhibits  low  Q  values  is  the  Tarim  block  in  central  Asia. 
The  low  values  there  might  occur  because  that  small  plate  is  located  within  the  Tethysides 
complex  or  because  the  resolution  in  our  inversion  is  not  sufficient  to  detect  it.  High  Q 
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values  in  the  South  China  Platform  are  not  unexpected  since  that  region  has  been  relatively 
stable  for  a  long  time. 

Conclusions 

Lg  coda  Q  throughout  Eurasia  exhibits  large  regional  variations,  between  about  1000  and 
200.  The  high  values  are  confined  to  northern  Eurasia,  in  a  long,  broad  band  that  extends 
from  western  Europe  through  northeastern  Siberia,  and  to  the  Indian  shield  and  South 
China  Platform.  Except  for  the  latter  two  regions,  all  of  southern  Eurasia  is  characterized 
by  low  Lg  coda  Q,  with  the  lowest  values  occurring  in  the  southernmost  regions.  These 
regions  include  those  that  have  been  most  intensely  deformed  by  the  Tethy sides  orogeny 
and  the  Arabian  Peninsula.  The  latter  region  was  intensively  deformed  in  its  eastern  portion 
while  its  western  portion  is  being  affected  by  process  occurring  during  the  spreading  of  the 
Red  Sea.  Low  Q's  in  the  Siberian  craton  lie  in  a  region  that  was  deformed  during  the 
Mesozoic  era. 

The  patterns  of  Q  variation  across  Eurasia  are  consistent  with  the  idea  that  Q  in  the 
upper  crust  is  directly  proportional  to  the  length  of  time  that  has  passed  since  the  last  major 
orogeny  in  any  region.  In  addition,  however,  the  observation  that  Q's  in  the  Germanotype 
Tethysides  are  higher  than  those  in  the  Alpinotype  Tethysides  suggest  that  the  intensity  of 
deformation  is  also  important  in  controlling  Q  in  the  upper  crust,  with  Q  in  more  intensely 
deformed  regions  being  lower  than  an  less  deformed  regions. 
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Figure  Captions 


1 .  Map  of  events  (+)  and  stations  used  for  inter-station  determinations  of 
phase  velocity,  group  velocity,  and  attenuation  in  the  Middle  East. 

2.  Single-station  surface  \vave  velocity  determinations  at  station  GNI  for  the 
earthquake  of  29  November  1994. 

3.  Single-station  surface  wave  velocity  determinations  at  station  ANTO  for 
the  earthquake  of  29  November  1994.\ 

4.  Inter-station  surface  wave  velocity  determinations  between  stations 
ANTO  and  GNI  for  the  earthquake  of  29  November  1994. 

5.  Inter-station  phase  and  group  velocities  for  the  station  pairs  in  Figure  1. 

6.  Inter-station  attenuation  coefficients  for  the  station  pairs  in  Figure  1 . 

7.  Top-  Synthetic  seismograms  for  fundamental  and  higher  modes  for  a 
crustal  model  of  the  eastern  United  States.  The  window  delineates  th  time 
interval  over  which  amplitude  spectra  were  determined.  Bottom  - 
Amplitude  spectra  computed  for  the  traces  above.  The  three  traces  and 
spectra  coirespond  to  vertical,  radial,  and  transverse  components. 

8.  Contoured  spectral  amplitudes  obtained  after  multiple  filter  analysis  of  the 
synthetic  seismograms  in  Figure  7.  The  largest  amplitudes  correspond  to 
the  expected  group  velocities  of  fundamental-mode  waves.  Smaller,  but 
still  easily  visible,  amplitudes  at  group  velocities  near  3.5  km/s 
correspond  to  higher-mode  surface  waves. 

9.  Polarization  parameters  (inclination  angle,  tilt  angle,  and  slope  angle)  and 
ellipticity  for  the  Rayleigh  wave  ground  motion  in  Figure  7. 


Introduction 

This  report  describes  our  work  on  fundamental-mode  surface  wave  propagation  in  the 
Middle  East.  For  that  work  we  are  proceding,  simultaneously  in  two  areas;  (1) 
Determination  of  phase  velocities,  group  velocities,  and  attenuation  coefficients  for  surface 
waves  in  the  Middle  East,  and  (2)  the  development  of  an  interactive  GUI  code  for 
determining  the  dispersion  and  polaration  of  seismic  waves.  Progress  in  each  of  these  areas 
is  described  in  the  following  paragraphs. 

Interstation  Determinations  of  phase  velocity,  group  velocity, 
and  attenuation 

We  searched  for  earthquakes  near  great-circle  paths  that  pass  between  various  pairs  of 
digital  broad-band  seismograph  stations  in  the  Middle  East.  Twelve  earthquakes  selected  in 
this  way  are  plotted  on  the  map  in  Figure  1.  Differences  between  the  directions  of  approach 
at  each  station  and  the  great-circle  path  between  stations  range  between  1.0  and  7.6 
degrees,  with  most  deviations  being  between  about  2  and  6  degrees.  Of  the  stations  shown 
in  Figure  1,  6  pairs  (GNI-ABKT,  ABKT-KIV,  GNI-ANTO,  ANTO-GNI,  ANTO-KIV, 
and  GNI-KIV)  have,  up  to  now,  yielded  useful  data.  Four  of  the  earthquakes  occurred 
beneath  the  Mediterranean  Sea;  for  that  reason,  paths  between  the  earthquakes  cross 
oceanic-continental  margins  that  may  laterally  refract  the  recorded  surface  waves  at  one  or 
both  stations,  causing  wave  propagation  directions  to  differ  from  the  deviations  given  in 
Table  1.  In  addition,  some  paths,  such  as  that  between  ANTO  and  KIV  and  between  GNI 
and  ABKT,  cross  regions  known  to  be  characterized  by  oceanic  crust  (in  the  former  case 
the  Black  Sea  and  in  the  latter  the  southern  Caspian  Sea). 

Examples  showing  the  determination  of  phase  and  group  velocities  for  Rayleigh  and 
Love  waves  between  stations  GNI  and  ANTO  appear  in  Figures  2-4.  This  is  a  path  through 
eastern  Turkey.  Figures  2  and  3  show  single-station  determinations  at  the  two  stations  and 
Figure  4  shows  the  interstation  results.  Similar  determinations  were  done  for  all  of  the 
events  in  Figure  1  and  Table  1. 

Figure  5  presents  12  interstation  determinations  of  Rayleigh  wave  phase  and  group 
velocity.  Large  variations,  especially  for  group  velocities  at  shorter  periods,  occur  from 
path  to  path.  At  longer  periods  velocities  are  similar  to  those  of  the  Basin  and  Range 
province  in  the  western  United  States  (Lin,  1989),  another  tectonically  active  region.  At 
periods  shorter  than  about  25  s  the  fastest  group  velocities  are  about  2.9  km/s,  values  that 
are  again  similar  to  those  in  the  Basin  and  Range.  Several  velocities  are,  however,  much 
lower  (as  small  as  2.2  km/s)  than  the  Basin  and  Range  velocities  of  Lin  (1989).  The  fastest 
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Figure  5 


Table  1.  Earthquakes  for  Two-station  Method 


Number  Date 

Origin 
h  m  s 

Lat 
(°  N) 

Lon 

(°E) 

Depth 

(km) 

mb 

Stations 

Az-dev  1 
(°) 

Az-dev  2 
(°) 

1 

24  May  94 

02  05  36.2 

38.66 

26.54 

17 

5.0 

GNI-ABKT 

7.1 

4.1 

2 

24  May  94 

02  18  34.9 

38.76 

26.60 

16 

5.0 

GNI-ABKT 

7.8 

3.9 

3 

03  Jun  94 

22  40  21.2 

28.74 

70.07 

33 

5.9 

ABKT_KIV 

6.9 

3.6 

4 

06  Oct  92 

08  57  17.7 

38.42 

56.52 

10 

4.7 

GNI-ANTO 

3.6 

1.9 

5 

21  Oct  94 

11  46  27.6 

38.25 

56.96 

33 

4.9 

GNI-ANTO 

4.1 

2.2 

6 

29  Nov  94 

14  30  28.4 

38.71 

20.48 

21  • 

4.9 

ANTO-GNI 

2.0 

1.0 

7 

26  Oct  93 

10  03  57.5 

36.74 

28.05 

33 

5.3 

ANTO-KIV 

6.3 

2.3 

8 

03  Jan  95 

22  51  45.0 

34.94 

23.55 

33 

4.9 

ANTO-KIV 

1.3 

9 

03  Feb  95 

22  29  13.3 

34.42 

25.03 

49 

4.8 

ANTO-KIV 

6.7 

3.3 

10 

30  Mar  95 

18  17  15.1 

34.40 

24.79 

7 

4.9 

ANTO-KIV 

5.8 

2.9 

11 

29  Mar  94 

07  56  53.9 

29.10 

51.26 

33 

5.4 

GNI-KIV 

7.6 

5.7 

12 

20  Sep  94 

05  51  46.0 

32.50 

48.77 

33 

5.0 

GNI-KIV 

4.1 

2.7 

13 

22  Apr  95 

00  21  48.6 

30.89 

49.91 

25 

5.1 

GNI-KIV 

5.8 

4.1 

36 


Rayleigh  wave  group  velocities  are  those  between  stations  GNI  and  KTV  west  of  the  Black 
Sea,  suggesting  that  the  crust  is  either  thinner  or  faster  than  normal  in  that  region. 

We  assume,  after  removing  effects  of  wavefront  spreading,  that  the  spectral  amplitudes 
(A)  of  seismic  surface  waves  attenuate  as  A(f)  ~  where  f  is  frequency,  y  is  the 
attenuation  coefficient  in  km'^,  and  x  is  distance  in  km.  The  attenuation  coefficients 
determined  as  a  function  of  period  appear  in  Figure  6.  As  with  velocities,  the  attenuation 
coefficient  values  exhibit  large  variations  from  path-to-path.  Mean  values  at  periods  of  30  s 
and  greater  are  about  0.4  x  lO"^  km  "k  At  shorter  periods,  they  increase  rapidly  with 
decreasing  periods  reaching  about  0.7  x  10"^  km'^  at  20  s  and  1.7  x  lO"^  km"^  at  10  s. 
These  can  be  compared  to  the  mean  interstation  values  of  Lin  (1989)  for  the  Basin  and 
Range  province,  a  region  of  crustal  extension  (with  moderate  tectonic  activity),  where 
seismic  wave  attenuation  in  the  upper  crust  is  known  to  be  high  (Lin,  1989;  Mitchell  and 
Xie,  1994).  Attenuation  coefficients  in  the  Basin  and  Range  are  about  0.3  x  10  km'^  at 
periods  greater  than  about  30  s,  0.35  x  10"3  km'l  at  20  s,  and  0.66  x  lO"^  km”l  at  10  s. 

Attenuation  coefficients  obtained  from  interstation  measurements  in  the  Middle  East  are, 
therefore,  higher  at  all  periods  than  those  determined  for  the  Basin  and  Range  using  the 
same  method.  Mitchell  (1995)  has  proposed  that  seismic  Q  (a  quantity  proportional  to  the 
inverse  of  attenuation)  in  any  region  is  proportional  to  the  length  of  time  that  has  passed 
since  the  last  episode  of  major  tectonic  activity  there.  For  both  the  Middle  East  and  Basin 
and  Range  to  conform  to  that  rule,  they  should  both  be  characterized  by  the  same 
attenuation  since  both  are  currently  active.  Since  attenuation  in  the  Middle  East  is  higher 
than  that  in  the  Basin  and  Range,  another  factor,  probably  the  greater  intensity  of 
deformation  in  the  Middle  East,  must  be  responsible. 

A  New  Tool  for  Surface  Wave  Dispersion  and  Polarization 
Analysis 

The  large  variations  of  surface  wave  velocities  and  attenuation  throughout  the  Middle 
East,  that  are  apparent  in  Figures  5  and  6,  suggest  that  our  measured  values  could  be 
significantly  biased  by  deviations  from  assumed  great-circle  paths  and  by  scattering  due  to 
lateral  complexities  in  the  crust  and  upper  mantle.  In  addition,  the  wave  forms  of  many  of 
the  surface  waves  that  we  found  sometimes  exhibit  choppiness  and  other  departures  from 
the  smooth  character  that  we  find  in  many  other  regions.  Because  of  these  complexities  it 
will  be  important  to  determine  the  particle  motion  of  the  surface  waves  that  we  study.  We 
have  recently  completed  coding  for  a  new  program  to  determine  surface  wave  dispersion 
and  polarization  as  a  function  of  period  or  frequency.  The  program  SPAN  (surface-wave 
polarization  analysis),  is  a  GUI  program  written  for  an  X-window  environment. 
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Figure 


Dispersion  can  be  determined  using  single-component  records  whereas  polarization 
measurements  require  records  of  three  components  of  motion.  The  computational 
procedures  follow  Levshin  et  al.  (1992),  but  our  implementation  differs  from  that  work  by 
permitting  researchers  to  proceed  interactively  through  the  process  using  results  of  each 
step  to  decide  on  parameters  to  be  used  in  the  next  step.  The  procedure  includes  (1)  plotting 
digital  time  series  and  displaying  a  dynamic  plot  of  particle  motion  within  any  selected  time 
window,  (2)  providing  a  view  of  the  selected  trace(s)  in  which  one  can  select  a  time 
window  of  interest  and  apply  selected  window  parameters  (length  of  time  window  and 
slopes  at  the  ends  of  the  window)  as  shown  in  Figure  7,  (3)  determining  the  Fourier 
spectra  for  the  selected  portions  of  all  traces,  also  shown  in  Figure  4,  (4)  displaying 
dispersion,  obtained  by  means  of  the  multiple-filter  method,  using  color  shading  for 
various  amplitude  ranges  (See  Figure  8  for  a  gray-scale  version),  and  (5)  plotting 
deviations  of  particle  motion  polarization  (including  inclination,  tilt,  and  slope  of  the 
polarization  ellipse),  again  with  color  shading,  from  that  expected  for  a  plane-layered  earth 
model  (Figure  8).  For  Rayleigh  waves  ellipticity  is  also  determined.  In  each  step  the  user 
controls  the  process  interactively  using  control  panels  in  the  active  window. 

Results  are  displayed  for  each  step  and  can  be  saved  and/or  printed  at  the  users 
discretion.  Messages  are  displayed  that  prompt  the  user  for  the  next  step  and  manual 
instructions  are  accessible  using  help  buttons.  Dispersion  curves  are  presented  as  contoured 
normalized  amplitudes  on  a  group  velocity-period  (or  frequency)  plot,  where  the  amplitude 
maxima  correspond  to  group  velocities.  Phase-equalization  and  time-variable  filtration  are 
used  to  isolate  the  desired  mode  for  improved  determinations  of  group  or  phase  velocity. 

Our  future  work  will  use  SPAN  to  routinely  determine  particle  motion  in  our  surface 
wave  analyses.  Surface  waves  that  are  found  to  deviate  greatly  from  a  great  circle  path  or 
exhibit  particle  motion  that  does  not  correspond  to  expected  Rayleigh  wave  or  Love  wave 
motion,  perhaps  because  of  mode  conversion,  will  not  be  used  in  our  modeling  of  velocity 
or  Q  structure  in  the  Middle  East. 
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